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Abstract 

The semiclassical formula for the quantum propagator in the coherent state representation 
{z"\e~^^'^^^\z') is not free from the problem of caustics. These are singular points along the 
complex classical trajectories specified by z', z" and T where the usual quadratic approximation 
fails, leading to divergences in the semiclassical formula. In this paper we derive third order ap- 
proximations for this propagator that remain finite in the vicinity of caustics. We use Maslov's 
method and the dual representation proposed in Phys. Rev. Lett. 95, 050405 (2005) to derive 
uniform, regular and transitional semiclassical approximations for coherent state propagator in 
systems with two degrees of freedom. 

PACS numbers: 02.30.Mv,03.65.Sq,31.15.Gy 
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I. INTRODUCTION 



Semiclassical methods are the fundamental tool in the study of the quantum-classical con- 
nection. In the limit where typical actions S become much larger than Planck's constant h, 
it is possible to use classical ingredients, usually classical trajectories, to produce approxima- 
tions to quantum mechanical objects, like matrix elements, wavefunctions and propagators. 
In Feynman's path integral approach to quantum mechanics, semiclassical approximations 
consist in realizing that the classical paths become dominant as S >> h and it suffices to 
add together the contributions of a small set of neighboring paths in the vicinity of the 
classical one. This apparently simple procedure, however, has two well known caveats that 
make the application of such formulas difficult: the existence of non-contributing classical 
solutions and the presence of focal points or caustics. 

The first of these issues, which is not going to be further discussed in this paper, is closely 
related to the Stokes Phenomenon, which is the abrupt change in the number of contributing 
solutions to an asymptotic formula when a certain boundary in parameter space is crossed 
Although a general criterion to decide whether a trajectory should be included 
or not as a true contribution to the formula exists, it is usually hard to verify in practice. 
An example of a careful study of these solutions can be found in j4|. More generally, one 
resorts to a simple a posteriori criterion: the contribution of each trajectory is computed 
and, if it leads to non-physical results, it is discarded. This kind of prescription have been 
widely used in the last years as, for example, in the semiclassical formula of the coherent 
state propagator in one [sl and two 6] spatial dimensions, in the momentum propagator [?] 
and in the semiclassical evolution of gaussian wave packets jsj. 

Singularities due to caustics is the other recurrent problem in semiclassical formulas. In 
the WKB theory ^] the semiclassical wave function in the position representation diverges at 
the turning points q = 0. In the momentum representation the equivalent problem occurs at 
the points where p = 0. In addition, for the Van Vleck propagator, which is a semiclassical 
formula of the propagator in the coordinate representation, (g"|e~*^^/^|g'), singularities 
occur at the focal points [l^. These are points along the trajectory from g(0) = q' to 
q{T) = q" where an initial set of trajectories issuing from the same initial point g(0) but 
with slightly different momenta, get together again, focusing at the same point q{t). 

The failure of the semiclassical approximation at these points, and the reason why a 
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singularity develops there, is that the usual quadratic approximation used to derive such 
formulas becomes degenerate and third order contributions around the stationary points 
become essential. The standard procedure to obtain improved formulas valid at caustics is 
due to Maslov 111 and it consists of changing to a dual representation where the singularity 



does not exist 
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12j. For a singularity in coordinates, one uses the momentum repre- 



sentation and vice- versa. The trick is that, when transforming back to the representation 
where the singularity exists, one should go beyond the quadratic approximation, otherwise 
the singularity re-appears. 

The subject of the present paper is the treatment of singularities due to caustic 
in the semiclassical formula of the coherent state propagator in two spatial dimensions 
K(z"*,z',T) = (z"|e~*''^^/^|z'). In spite of the fact that this is a phase space representation, 
where no turning points exist, this propagator is not free from caustics j^, 0, [l3, 14 1, al- 



though earlier works on the subject indicated so These points have been 

termed phase space caustics. 

The caustics in Ksc(z"*, z', T) have the same origin as the focal point divergence in the 
Van-Vleck propagator, namely, the breakdown of the quadratic approximation. Therefore, 
it is natural to seek for a dual representation as in Maslov's method to derive higher order 
approximations. However, since both coordinates and momenta are used in the coherent 
states, there seems to be no room for a natural dual representation. In a recent paper |l9| 
we have proposed the construction of an application between f{z*) = {z\tp) and an associate 
function f{w) that plays the role of the dual representation for the coherent state propa- 
gator. Using this auxiliary mapping we were able to derived a uniform approximation for 
the propagator of one- dimensional systems that is finite at phase space caustics. In this 
paper, we use such a representation to derive regular, transitional and uniform semiclassi- 
cal approximation for the coherent state propagator of two-dimensional systems, which is 
the simplest case where conservative chaos is possible. The resulting formulas involve, as 
expected, the Airy function and the third derivatives of the action function. 

This article is organized as follows: in Sect. II we review the semiclassical formula for the 
coherent state propagator in two dimensions and discuss its singularities. In Sect. Ill we 
review the dual representation proposed in Ref. 19[] and generalize it for two-dimensional 
systems. In Sect. IV we use this representation and the Maslov method to derive regular, 
transitional and uniform approximations for the coherent state propagator. Our conclusions 



3 



and final remarks are presented in Sect. V. 



II. THE SEMICLASSICAL LIMIT OF THE COHERENT STATE PROPAGATOR 

In tliis section we briefly discuss the usual semiclassical formula for the propagator in 
the coherent state representation. The 2-D non-normalized coherent state |z) is the direct 
product of two 1-D states, |z) = ® \zy), where 

\z^) = e^'-"^|0), 

1 / .pr 

Zr = —;={—+ I — 
a/2 \br Cr 

The index r assumes the values x or y. |0) is the ground state of a harmonic oscillator 
of frequency Ur = Ti/{mh1), a\ is the creation operator and g,., pr are the mean values of 
the position q^. and momentum p^- operators, respectively. The widths in position h^. and 
momentum Cr satisfy hrCr = h. In addition, the complex number Zr is eigenvalue of with 
eigenvector \zr). 

The coherent state propagator K(z"*,z',T) = (z"|e~*^"^/^|z') represents the probability 
amplitude that the initial coherent state |z') evolves into another coherent state |z") after 
a time T, according to the Hamiltonian H. Notice that, since the initial and flnal coherent 
states are non-normalized, all the propagators considered in this paper should be multiplied 
by e~^'^''~5l^"l^ to get the usual propagators with normalized bras and kets. 

The semiclassical approximation for this propagator was flrstly considered by Klauder 
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221 ] and Weissman 23|] . More recently, however, a detailed derivation was presented 



or systems with one degree of freedom 



24| . The expression for two-dimensional systems is 



4?(."-.z'.T) = 5:.U_exp 1^^, (2) 



traj. 



IdetM^^I h 



where the index ^"^^ was inserted to indicate explicitly that this formula was obtained by 
means of a second order saddle point approximation. The factors Mw and depend on 
(generally complex) classical trajectories. These trajectories are best represented in terms 



of new variables u and v, instead of the canonical variables q and p, defined by 

1 / ?r , .Pr\ , 1 / .Pr 



The sum in Eq. ([2]) runs over all trajectories governed by the complex Hamiltonian H{u, v) = 
(v|i7|u). They must satisfy the boundary conditions u(0) = u' = z' and v(T) = v" = z"*. 
Notice that Qr and Pr are complex variables, while the propagator labels {q[., p'^ for the initial 
state and g", p" for the final one) are real. In Eq. ([2]), JF is given by 

^(v", u', T) = 5(v", u', T) + GiV, u', T) - l^v., (4) 

where S, the complex action of the trajectory, and Q are 

iS(v", u', T) = / — (u V — u v) — c^t — — [u"v" + uV] , (5) 




The matrix Mw is a block of the tangent matrix defined by 

(7) 



the tangent 
i). Finally, 





where 5u and Sv are small displacements around the complex trajectory. We use a single 
(double) prime to indicate initial time t = (final time t = T). The elements o: 
matrix can be written in terms of second derivatives of the action (see Ref. 
(Jvv is the phase of det Mw 

The classical trajectories contributing to the propagator are functions of nine real pa- 
rameters: four initial labels q'^, q'y, p'^ and p'y, four final labels q", q'y, p" and py, and the 
propagation time T. As one changes one of these parameters continuously, it might happen 
that two independent solutions become very similar to each other. In the limit situation they 
might coalesce into a single trajectory, characterizing a phase space caustic, or a bifurcation 
point. At the immediate neighborhood of the caustic these solutions will satisfy identical 
boundary conditions. Therefore, close to the caustic, we can set small initial displacements 
6u' = and 6V 7^ in such a manner that, after the time T, 6u" ^ and 6V = 0. Eq. ([7]) 
then reduces to 



6u" \ IM^^M^^\ 
) \ M,, M,, ) \ 6V 



(8) 



implying that det Mw = 0. Therefore, at the caustic the pre-f actor of the Eq. ([2]), 
I det Mvvl"^''^, diverges and the semiclassical formula cannot be used. The main purpose 
of this paper is to correct the semiclassical formula in these situations, replacing the singu- 
lar pre-factor by a well behaved Airy-type function. 

As mentioned in the Introduction, in addition to the divergence of the pre-factor, the 
semiclassical formula ([2]) can exhibit other peculiarities, which we shall not address here. For 
example, for some complex trajectories the imaginary part of JF can be large and negative, 
giving unreasonably large contributions to the propagator. This problem is related to the 
Stokes Phenomenon, and lead to the exclusion of some trajectories from the sum in Eq. ([2]) 
in order to eliminate the unphysical results they produce 



III. DUAL REPRESENTATION FOR THE COHERENT STATE PROPAGATOR 

The main difficulty in dealing with phase space caustics is the lack of a dual representation 
for the coherent state propagator. Caustics in position representation can be removed by 
changing to the momentum representation and vice- versa. For the coherent state propagator, 
since both position and momentum are being used, it is not clear what to do to bypass the 
phase space caustics. In Ref. [l^ we defined an application that plays the role of a dual 
representation for the coherent state basis and used it to find a uniform formula for the 
semiclassical propagator for one-dimensional systems. The purpose of this section is to 
extend the formalism of Ref. for systems with two degrees of freedom. 

Based on the relations 

which can be demonstrated by differentiating the complex action of Eq. ([5]), we per- 
form a Legendre transform on 5 (v", u',T) replacing the independent variable v" by 
u" = {i/h){dS /dV). The transformed function S depends on the variables u' and u", 
instead of u' and v", 

S (u", u', T) = S (v", u', T) + ihu'V , (10) 



and satisfies the relations 



„ i dS , i dS 

h ou" n ou' 
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With these properties in mind we define a dual representation K (u", u',T) for the propa- 
gator K (v", u', T) by 

K(u", u',T) = — [ K(v", u',r)e-""""dV', (12) 



27r 



c 



K (v", u',T) = ^ I K (u", u', T) e^"^"d'u", (13) 
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where C and C are convenient paths that, as specified in pL3], are chosen in such a way that 
Eqs. f|T2|) and f|T3|) are a Laplace and a Mellin transform, respectively. The analogy between 
these two expressions and the corresponding coordinate and momentum representations is 
not complete. This is because, while K(v", u', T) is the quantum propagator, K(u", u', T) 
does not seem to correspond to an explicit quantum matrix element. It is a mathematical 
artifice that allows for the continuation of the propagator in an auxiliary phase space, rather 
than a quantity with a direct physical interpretation. 

In order to obtain a semiclassical formula for K (u", u',T), we insert Eq. ([2]) into (|T2l) . 

Ksc(u", U',T) = — / e^'5(^"."''^)+ie(^"."''^)-Wv-|ln|detM..|-u"v"^2^//_ ^^^^ 

27r Jc 

Rigorously, Eq. (JT4l) says that to calculate Ksc for a set of parameters u", u' and T, we 
need to calculate the contribution of the trajectory beginning at u(0) = u' and ending at 
v(T) = v", and sum over all v" lying in the path C. Notice that, for each trajectory, the 
value of the variable u at time T is function of u', v" and T, namely, u(T) = u(v", u',T). 
In the semiclassical limit this integral can be solved by the steepest descent method (3|, 
according to which the critical value v" of the integration variable satisfies 

I dS 



or u" ^ 



(15) 



where we have considered that Q and In | det Mwl varies slowly in comparison with S, since 



the former is of order h while the later is of order (see Ref. 2J]). Eq. dlSI) says that the 
critical trajectory satisfies u(0) = u' and u(T) = u(v^',u',T) = u", i.e., the critical value 
v^' of the integration variable is equal to v(T) of a trajectory satisfying these boundary 
conditions. This shows that the integration path C must coincide with (or be deformable 
into) a steepest descent path through v". Expanding the exponent up to second order 
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around this trajectory and performing the resulting Gaussian integral we obtain 



K«(u".u'.T)=5: J^gjL^ exp{l5(u''.u'.T) + ia(u''.u'.T)-l.„.} . (16) 

traj. y I I ' 

where, again, the index ^"^^ is used to indicate the method of integration used. The sum over 
stationary trajectories was included because more than one of them may exist. To derive 
the last equation, we have also used the result 

- det I ^''^ ] = ;.^[^^^j^-|e-(— ), (17) 



'-VV 



with Sai3 = d^S/dadp, for a,P = v" or Vy, and (Xuv is the phase of detMuv This last 
equation can be obtained by considering small variations of Eq. ([9]), rearranging the terms 
so as to write Su" and 5v" as function of 5u' and 6v', and comparing with Eq. ([7j). 

The new semiclassical propagator Ksc is a function of complex classical trajectories sat- 
isfying u' = u(0) and u" = u(T). Muv is given by Eq. (JTj), Q (u", u',T) is the function Q 
calculated at the new trajectory, and S {u" , u',T) is given by Eq. ( ITOl) . It is easy to see 
from Eq. ([8]) that, when detMw is zero, detMuv is not. This is a fundamental property 
that one has to bear in mind when deriving approximations for K (v", u',T) inserting Kg^ 
into Eq. f|T3|) . Three such approximations will be derived in the next section. 



IV. COHERENT STATE PROPAGATOR FROM ITS DUAL REPRESENTATION 



Replacing Eq. (fT6|) back into Eq. (fT3|) we obtain 



Ksc(v", u',T) = — / e^'5'(""'"''^)+^^(""'"''^)-5'^"— ii>^l<i'=*^^"-l+""^" d^u". (18) 
27r Jc 

To solve Ksc for the parameters v", u' and T, we need to sum the contributions of all 
trajectories beginning at u' and ending at u" lying in C. The saddle point u" of the 
exponent satisfies 



d 



o — tnu V 



or 



(19) 



which says that the most contributing trajectories are those with boundary conditions 
v(T) = v" and u(0) = u', exactly as in Eq. ([2]). Therefore, expanding the exponent up 
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to second order around the critical trajectory, solving the remaining Gaussian integral, and 
using the result (see Eq. flA22l) of the appendix) 

- det ( ]^-det S„... = /.^If ^-|e-(-— ), (20) 

we recover Eq. ([2]). 

Clearly, the connection between the propagators of Eqs. ([2]) and (flBl) via steepest descent 
approximation with quadratic expansion of the exponent works only in the regions where 
both detMuv and det Mw are non-zero. Close to caustics, where det Mw = 0, Kg(^ is still 
well defined and Ksc can be obtained by doing the inverse transform f|T8l) but expanding 
the exponent to at least third order. There are, however, several ways to handle such an 
expansion, depending on how close to the caustic a given stationary trajectory is. In the 
next subsections we show how to obtain three approximate formulas for the propagator: 

In Sect. IIV A[ we evaluate Eq. flTHl) by expanding its integrand up to third order around 
the stationary trajectories. As a result we find that each contribution already present in 
Kgc appears multiplied by a correction term X^. This regular formula for the semiclassical 
propagator is good only if the stationary trajectories are not too close to caustics, so that 
second and third order terms contribute to the integral. 

In Sec. IIVBI we consider the situation where two contributing solutions are so close each 
other that, if we used the regular formula, the contributions would be counted twice. We 
therefore perform a transitional approximation, where the exponent of (ITS!) is expanded 
around the trajectory that lies exactly at the phase space caustic. Since this trajectory is 
not generally stationary, this approach works only if the stationary solutions are sufficiently 
close to the caustic. 

Finally, in Sect. IIV CI we derive a uniform approximation, which is applicable both near 
and far from the caustics but might not be so accurate as the two previous expressions. 



A. Regular Formula 

The philosophy of the regular approximation is to correct the contribution of each station- 
ary trajectory by including third order terms in the expansion of the exponent of Eq. (ITSl) . 
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When this expansion is performed we obtain 
k(^^ (v", u', T) = 



det M. 



e-k^ X In (v", u', T) 



(21) 



where the quantities between brackets are the same as in Eq. (j2]), and the correction term 
Xr is given by 



1r 
and 



det Su"u" 



D = E = ^Su'^u'^u'^, F = \Su'^u'j,u'j. and G = ^Sui',u',',u'j,- 



(23) 



6 ""y ""y ""y 

All functions and constants in Eq. ( 12T1) are calculated at the critical trajectory. In Eq. (1231) . 
we define Sa/s^ = {d^S / dad[3d'^) and S^p = (d'^S/dadP), for 7 = or Uy. The 



integration contour of Eq. fl22|) is chosen to coincide with the steepest descent of the saddle 
point. 

The integral fl22|) has no direct solution. However, it can be largely simplified in the 
coordinate system that diagonalizes the matrix the quadratic terms, 



A B/2 
B/2 C 

Therefore, we perform the change of variables 



2Su"u" 



(24) 





1 




-iV_(A- A+) 




5uj^ 




i?/2(A_ - A+) 


N+B/2 


-N_B/2 




6u_ 



(25) 



where N± are normalization constants and \± are eigenvalues of |Su"u", 



N± = ^ (5/2)' + {A- X±y and A± 
In the new variables Eq. (1221) becomes 



1 ± 



1-4 



det S, 



\ (trSu"u"y 



(26) 



R 



A+A_ 



j <\[5u+](l[5U-] (.Ti{^+^^l+^-^^''-+D'5ul+E'5ul&u.+F'5u+&ul+G'5ul} 



(27) 
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where the new coefficients, D' ^ E', F' and G", are combinations of those in Eq. (123|) . Our 
final formula depends just on G", which amounts to 



G' 



N. 



A- Ah 
B/2 



D 



A- Ah 
B/2 



E + 



A- Ah 
5/2 



F + G 



(28) 



According to Eqs. (!20|) and (!24l) . when det Mw ^ 0, det Su"u" also tends to zero, causing 
the breaking down of the quadratic approximation. However, in terms of the variables 6u+ 
and we see that det Su"u" (= 4A+A_) goes to zero in a particular way: while A_ vanishes. 
Ah- generally remains finite. Therefore Eq. (1271) is always of a gaussian type integral in the 
variable Su+, since we are still able to neglect third order terms in this direction. Solving 
the integral in 6u+ leads to 



-zA_ 



(29) 



3G' \ V3 



Now we perform a last changing of variables, t = (^) 



and obtain 



2v^ 



(30) 



where w = and U{w) is given by 



27T 



wt + 



(31) 



for i = 1,2,3. The index i refers to three possible paths of integration Cj, giving rise to 
three different Airy's functions (see Ref. [s!])- Rigorously the choice of the path should be 
done according to Cauchy's Theorem, after all the path to be used has to be obtained by 
a deformation of the original contour of integration. In practice, however, it might be very 
difficult to find the correct path in this way, and we have to use physical criteria to justify 
the choice of Cj. 

Inserting (l30l) into Eq. ( 1211) and considering the existence of more than one critical tra- 
jectory, we finally find the regular formula 



traj. 



det M. 



g^^(v",u',T) 



2y/n w e-^"" fi(u;' 



(32) 



In this equation, each stationary trajectory gives a contribution which is that of the 
quadratic approximation multiplied by a correction factor Jr that depends only on the 
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parameter w. Close to a caustic A_ is very small but G' (generally) remains finite. Exactly 
at the caustic \w\ is zero, getting larger and larger as we move away from it. Therefore we 
expect that should go to 1 as \w\ goes to infinity, since the regular expression should 
recover Kg^; in this limit. To verify this assertion, we look at the asymptotic formulas for 



the Airy's functions 25 1 



fl 


» 




l/4g- 


2^3/2 
3 


k 






l/4g| 


U,3/2 


h 


» 


-27^^" 


l/4g2 





(33) 



Using these expressions in Eq. ( [30l) . we see that only hiw) produces the desired asymptotic 
result, indicating that this is the proper choice of Airy function far from the caustic. However, 
this is so only because we have taken the principal root in the definition of w. As the physical 
results should not depend on the arbitrariness of branches in the complex plane, the choice 
of a different root would lead to a different path Cj, so that physical results remain the same. 
A careful discussion about this point can be found in {4]. 
Exactly at the caustic, w = 0, the regular formula becomes 

where the sum was excluded because the critical trajectories coalesce at this point. 



B. Transitional Formula 



Each contribution to the semiclassical propagator calculated in the last section (as well 
as those of Eq. ([2])) has information about the critical trajectory plus its vicinity. If two 
trajectories are very close each other, like in the vicinity of a phase space caustic, their regions 
of influence might overlap. The regular formula cannot be used in these situations, since it 
assumes that the trajectories can still be counted independently. To find an approximation 
for K (v", u', T) valid in this scenario, we shall perform the integral f|T8|) expanding the 
exponent about the (non-stationary) trajectory corresponding to the phase space caustic 
itself, defined by the point u" where | detSu"u"| and therefore | detMwl [see Eq. fl20|) ] are 
zero. Evaluating this single contribution to third order should be equivalent to include and 
sum over each stationary trajectory. 
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The condition det Su"u" = leads to the PSC trajectory that begins at u(0) = u' and 
v(0) = v', and ends at u(T) = u" and v(T) = v", where v" is assumed to be close to v". 
Expanding the exponent of Eq. (JTSl) up to third order around this new trajectory yields 



K 



PSC 

sc 



(v", u', T) = (detM„0"'^'ei^^(^"'"'^)+^(^"'"'^)-^'^^"(^"-^")> Ij 



(35) 



where 



27r 



with X = dS/du" —ihv'^ and Y = dS/dUy —ihVy. The functions appearing in Eq. (13^ and 
all the coefficients are calculated at the PSC trajectory. 

We solve Eq. fl36l) using the same technique described in the last section, with the use 
of the transformation (!25|) . However, as we deal with the PSC trajectory, A+ = A + C and 
A_ = 0. The integral Xj- becomes 



Xt 



2lT 



(37) 



where the only coefficients that appear in the final formula are 



A - A_ 
B/2 



X + Y 



A+-A_ 



A- Ah 
B/2 



X + Y 



(38) 



and G', given by Eq. (128|) . 

The integral over 6uj^ can be performed neglecting terms of third order. We obtain 



Xt = e~^^ f d[5u_]exp|^ [b6u^ + G'5ut] \. 



(39) 



By setting t = (^)^^^ 6u_, the last equation can be written as 



Xt 



'ink 

a7 



3G" 



(40) 



where w = (ggf^^-ji/s and the function fi{w) refers to the Airy's functions (13T|) . Finally, we 
write the transitional formula by combining Eq. fHUl) with Eq. (15^ . 



ih 



TX 



1/3 



_i_ a 



fi(w) e 



|[5+e]-i?iu"(v"-v") 



(41) 



A+(detMuv) VSG"^ 

Equation (HT]) depends on the PSC trajectory, which satisfies u(0) = u' and v(r) = v", and 
is valid only if v" is close to v". 
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Far from the caustic Eq. (HT!) does not make sense, since the PSC trajectory becomes 
completely different from the actual stationary trajectories. On the other hand, when the 
propagator is calculated exactly at the PSC, Eqs. (l^Tj) and ([22]) should furnish the same 
result. This can be verified by setting v" = v" and a = b = w = Oin Eq. fHT]) . which 
reduces directly to Eq. (^^. 



C. Uniform Formula 



The regular formula is good as long as one is not too close to a phase space caustic, 
whereas the transitional formula is good only very close to it. In either cases the expressions 
we derived cannot be used everywhere in the space spanned by the parameters u', v" and 
T. The uniform approximation provides such a global formula [20]. The basic idea is to 
map the argument of the exponential in (fTS!) into a function having the same structure of 
saddle points as the original one, i.e., two saddle points that may coalesce on the phase 
space caustic depending on a given parameter. 

In order to simplify our calculation, we shall use the variables and u'!_, instead of the 
original u'^ and Uy [see Eq. fl25|) ]. In these variables the exponent of Eq. f|T8l) 

E (u", u', T) = (u", u', T) + (u", u', T) - '-a^^ - ^ In | det M.^l + u"v", (42) 

becomes 

S{ul,u'^)=E [u" {ul , m': ) , u' , T] , (43) 

where we omit the dependence on the variables u' and T because they are not being inte- 
grated. The integral flTSl) then becomes 



1 

2^ 



e^K'"-^ du'[du'^. (44) 



Since the main contributions to this integral comes from the neighborhood of the saddle 
points, we can map the exponent £{u'\_,u"_) into a new function N{x,y), where x = x{u'\_) 
and y = y{u"S). We restrict ourselves to the case where there are only two critical points, 
u" = (u",u")i and U2 = (M",n")2, which, depending on the parameters u' and T, may 
coalesce at the phase space caustic. Then 

L / e'^^'""-) du"Au"_ ^ — [ J(x,i/,)e^("'^) dxdy. (45) 
2n J 2n J 
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The simplest function with these properties is 



N{x, y) = A- By + — + Cx^ 



(46) 



where A, B and C may depend on u' and T. The mapping requires that the saddle points 
of N{x,y), which are (0, ±VB), coincide with the critical points u'l^'- 



S{u'[) =£i = Ar(0, ^/B)=A- lB^'\ 
£(u'2') =82 = iV(0, -VB) =A+ lB^I\ 



implying that 



A=]^{£i + £2) and B 



2/3 



(47) 



(48) 



Another condition required to validate the method is to impose the equivalence between the 
vicinity of critical points of N{x,y) and £{u'^,u"_), 



5N + + -d'^N + 
2 6 



5£ 



2 6 



(49) 



{o,±Vb) 

This equation allows us to find how to transform an arbitrary infinitesimal vector (5'u" , ) 
into {6x,6y), around the critical points. It provides, therefore, information about the Jaco- 
bian J{x, y) of the transformation calculated at the critical points, namely, Ji = J(0, \/B) 
and J2 = J{0,-y/B). 

As the first derivatives of S and vanish at the critical points, Eq. fj49l) implies that 



-[ 6x 6y 



dxdy dydx"^ 
dydx dxdy^ " dy^ 3 dy^ " 



{0,±v^) 




(50) 



should be equal to 



2 ( 5u'L 



du" du'i 



1 d'^^e 



3 du- 



5u'L 



- + 



+ 



au'idu"_'^ 



5u" 



d^8 

1 d^e 



6u'' 



) 








1 





(51) 



Writing (J-u'l = a+6x and ^m" = a_6y this equality results in 



2C, 



I a^£ 



a^e 



au" au'' 



+ 



du'tau" 



0, 



■■1,2 



(a_5?/) 



au' 



1 



^ + 39^?^ 



[a_6y) 



a+a_ j> = 0, 

±2v^ + Uy. 



(52) 
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In the limit of small h, Q and det Muv vary slowly in comparison with S and the first and 
last of equations fl52l) become, respectively, 



(53) 



l{[K + D'{a+5x)]a\]\^,,^=C, 
t { [A_ + G' {a^5y)] al]\^,,^ = ±VB+\5y. 

Moreover, the second and third fl52l) imply that E' = F' = 0. We emphasize that D', E', F' 
and G' are the same coefficients as those of Sec. IIV A[ 

Eqs. ( I53i) can be solved if we neglect the terms containing Sx and Sy. We find 



-ihC 



and 



so that the Jacobian at the saddle points becomes 



(«-)lu^ 



(54) 



Ji,2 = ia+a_)\^,. 



(55) 



The full Jacobian can therefore be conveniently written in the vicinity of the saddle points 



as 



and the uniform approximation for the propagator becomes 

KUN(v",u',T) = 1-1 Jix,y) e^-'^y^y'/'^'^'dxdy. 
Performing the integral over x we obtain the final expression 



Ksc (v ,u ,T) = i^-n > ' 
where fi is given by Eq. (|3T|) and 



B 



(56) 



(57) 



(58) 



\ (4A+A„)|„,,^ 



Vb 



det 

det M. 



(59) 



Eq. ( l58l) is the uniform formula for the two-dimensional coherent state propagator. As 
in Sects. IIV Al and IIVBI the determination of the proper path of integration Gi is done by 
physical criteria. 

Eq. (!59|) shows us how the singularity in the coalescence point is controlled. When 
det Mvv goes to zero, the difference between Si and S2 also vanishes, so that the quotient 
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\fBI det Mvv [see also Eq. (HSIl ] remains finite. Notice, however, that this fraction might 
become extremely fragile close to a caustic, because both numerator and denominator go 
to zero. Exactly at the caustic we can return to the second of Eqs. fl53l) to find the correct 
value of a_: 

One should also remember that, if h is not sufficiently small, the derivatives of Q and det Mw 
may become important, specially when A_ 0. 

It is interesting to check that the uniform approximation ( l58l) recovers the quadratic 
approximation away from the caustics, i.e., in the limit B oo. According to Eqs. ( l33l) we 
find that, for large w, 

(61) 



w~ 


-l/2f/ 


[w) 




_l/4g-|^3/2 


w~ 




[w) 






w~ 




» 


~ 27^^" 


-l/4g§«;3/2^ 



Inserting Eqs. fl55]) and flM]) into the uniform approximation results in 

' 2 2 

—ig2e-^~i^ B~^^^, by using fi 

K^c^(v",u',T)^ 



-(7le^+i^'^'i3-l/^ by using fa • (62) 

gie-^+^'^'^^^B~'^^^, by using fs 



It's easy to see that using the contour Ci + C2 we find |Kg^^| = |Kg^|. 

Another way to arrive at the same conclusion is as follows: if u'/ and are not close 
each other, we can individually evaluate the contribution of each one through the second 
order saddle point method and sum the contributions at the end. Starting from Eq. (!57|) we 
get 

KUN(v",u',r) = ^JJ{y)e^-^y+y'/'dy 

Finally we consider the uniform formula evaluated exactly at the caustic. To do so we 
rewrite Eq. fl571) using the uniform Jacobian given by Eq. fl60l) : 



KUN(v",u',T) = i- 
Zn 



Since ^-^^ J ey'/^dy = e-^-^/^ fi(0), we find the same result as found previously with the 
formulas of the Sects. IIV Al and IIVBI calculated at phase space caustics. 
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V. FINAL REMARKS 



Semiclassical approximations for the evolution operator seem to be plagued by focal points 
and caustics in any representation. A relatively simple way to derive improved expressions 
that avoid the singularities of such quadratic approximations is provided by the Maslov 
method. The method explores the fact that, for example, the coordinate representation of 
the propagator, {x\K{T)\x') can be written as the Fourier transform of the propagator in 
its dual representation, {x\K(T)\x') = f {x\p){p\K(T)\x')dp. If the trajectory from x' to 
X in the time T corresponds to a focal point, we can still use this integral expression and 
the usual quadratic approximation for {p\K{T)\x') , as long as we perform the integral over 
p expanding the exponents to third order around the stationary point. This results in a 
well behaved approximation for the coordinate propagator in terms of an Airy function. In 
this paper we have shown that a similar procedure can be applied to the coherent state 
representation and derived three similar third order formulas that can be used depending 
on how far the stationary trajectory is from the phase space caustics. Although we have 
considered only systems with two degrees of freedom the extension to higher dimensions is 
immediate. We note that a uniform formula for the coherent state propagator was previously 
derived in 27| for a particular Hamiltonian. 



The regular formula fl32l) is the simplest of our three approximations and consists of 
a sum over the same complex trajectories that enter in the quadratic approximation. The 
contribution of each trajectory is regularized by a term that avoids divergences at phase space 
caustics. We emphasize that this regularization deals just with the problem of caustics, so 
that we still need to identify contributing and non-contributing trajectories in order to get 
acceptable results. This approximation holds as far as the contributing trajectories are not 
too close to the caustics, otherwise the vicinities of different trajectories can start to overlap 
and their contributions would be miscounted. The transitional formula fj^Tj) works exactly 
in this situation. It involves the contribution of the PSC trajectory alone, and therefore is 
valid only very close to the caustics. Finally, the uniform formula ( l58l) is valid everywhere, 
near of far a caustic. The formula we derived deals with the simplest topology of caustics 



28|. 



All three semiclassical formulas derived here involve the calculation of third order deriva- 
tives of the action. We presented an algorithm to evaluate these derivatives numerically in 
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Appendix |X1 Numerical results using these expressions will presented in a future publication. 



APPENDIX A: DERIVATIVES OF THE ACTION S 

In this appendix, we show how second and third derivatives of 5(u',u",T) can be cal- 
culated for a given trajectory. This procedure can be used with any set of variables (for 
example, (u', v",T) or (q',q",T)) with minor modifications. 

1. The Tangent Matrix and the Tangent Tensor 



The equations of motion in the u and v variables can be written in compact form as 



f. = 7 H' 



where the vector r and the matrix J are given by 

( u,,\ / -i/n \ 

and J = 



(Al) 



and 



H 



-i/n 

i/h 
\ i/h J 



dH 



(A2) 



(A3) 



Expanding Eq. flAl|) up to second order around a reference trajectory r(t), we obtain 



Sr-i = Jij H-^ 5rk + -Jij 6ri Hi^^ 5r, 



where 



and H, 



Ikj ' k , 



(A4) 



(A5) 



The solution of Eq. flA4|) can be expressed in terms of the initial displacement 6r{0) as 



6r,{t) = M.,,{t) 5r,(0) +5r,(0) Uku{t) 6n{Q) 



(A6) 
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where the tangent matrix M and the tangent tensor U satisfy M(0) = 1 and t/(0) = 
0. Differentiating this equation with respect to t and by using Eq. flA4p . we obtain the 
differential equations satisfied by M and U directly: 

\ ) 

where we have discarded terms of third order in (5rj(0). This leads to 

= Ju Hi Mkj (A8) 

and 

Uijk = Jkl Uijm + -Jki Mni H^^i Mmj ■ (A9) 

These two sets of differential equations can be solved for a given reference trajectory f{t) 
and boundary conditions M(0) = 1 and U{0) = 0. 



2. Derivatives of S 



Here we show how to obtain the second and third derivatives of S in terms of M and U. 
We start from Eqs. (ITT]) , which can be written as 



(AlO) 



where 

< 



V 



f i/n \ 

i/h 
-i/h 
V -i/h ) 



(All) 



and S\ = dS/dUi. 

Considering variations on Eq. flAlOp around the reference trajectory and expanding up 
to second order, we get 

1 

2' 

where 



5V, = K,, S^, 5Uk + -K,, 5Ui S^,^ 6Uk , 



(A12) 



s: 



and S, 



(A13) 
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The idea now is to manipulate Eq. (1A12|) so that final displacements are written in terms of 
the initial ones. To do this we write 



6V = A6r(0) + B6r(T) 

^ ^ ^ ^ (A14) 

SV = C5r(0) + D6r{T) 

where A, B, C and D are 4x4 matrices that can be written in terms of 2 x 2 blocks as 

"^(ooj' "^(loj' "^(ooj' ^^(oij 

Replacing Eqs. flA14p into (lA12p and solving for 6r(T) produces 

6r{T) =[d-K S" [k S" A-C^ 5r(0) + ^A-^ w, (A16) 

where A = K-\D - KS"B) and 



= [A6r{0) + B6r{T)], S^, [A6r{0) + B6r{T)U 
^ [A6r{0) + BM5v{0)]i S^^ [A5r{0) + BM5r{0)], 
= [L6rmi S'L [L6vmm ■ 



(A17) 



In this expression we have discarded terms of third order in 5r(0) and we have defined the 
auxiliary matrix L = A + BM. Computing all these matrices explicitly, we find 

D-KS"B] '= ( ° I, (A18) 

'-'u"u" '-'u'u" ^ 

KS"A-C)=\ I , (A19) 



L-'=\ ^ " I , A = ° I ■ (A20) 



Comparing linear terms of Eq. (]A16[) with ( ]A6l) . we find 



M = I = ( . , _ I (A21) 

Mvu Mvv / \ (V^) ( 'S'u"u" S^i^u ^u'u' — 5'u"u' j — 'S'u"u" S^,^„ 
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or 



Finally, comparing the quadratic terms, 



= Ukii Sri 



(A23) 



or 

-Aj/ Lnk S^mj Lml = Ukli ■ (A24) 

Solving for the third derivatives of 5* produces 

^ijk = 21/^1 Afcn Umln L^j^, (A25) 

where A and are given by ( ]A20I) . 
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